Phase-Sensitive Detection for Unconventional Bose-Einstein Condensations 
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We propose a phase-sensitive detection scheme to identify the unconventional px ± ipy symme- 
try of the condensate wavefunctions of bosons, which have already been proposed and realized in 
high bands in optical lattices. Using the impulsive Raman operation combining with time-of-flight 
imaging, the off-diagonal correlation functions in momentum space give rise to the relative phase 
information between different components of condensate wavefunctions. This scheme is robust 
against the interaction and interband effects, and provides smoking gun evidence for unconventional 
Bose-Einstein condensations with nontrivial condensation symmetries. 
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Macroscopic condensates of bosons and paired 
fermions are of central interest in condensed matter 
physics. Order parameters of Cooper pairings are termed 
"unconventional" if they belong to non-trivial represen- 
tations of rotational symmetry groups. Celebrated exam- 

)les include the d-wave pairing states of high Tc cuprates 
the p-wave pairing states of ^He A and i?-phases 

i, and Sr2Ru04 0, Among various experimen- 
tal tools, phase-sensitive detections are exceptional as 
they provide the smoking-gun evidence for unconven- 
tional pairing symmetries, such as the 7r-phase shifts in 
the joint corner SQUID junctions [l| and the tri-crystal 
superconducting ring experiments [4I of high Tc cuprates. 
Unconventional symmetries have also been generalized to 
the particle-hole channel pairings, i.e., the Pomeranchuk 
type Fermi surface instabilities in both density [6|] and 
spin channels 7, 8] . The spin instabilities in high orbital 
angular momentum channels are denoted as "unconven- 
tional magnetism" . 

Recently, unconventional symmetries have been intro- 
duced to the single-boson condensates , denoted as 
"unconventional" Bose-Einstein condensations (UBEC). 
Their condensate wavefunctions belong to non-trivial 
representations of the lattice point group. Such states 
have been proposed in high orbital bands of optical lat- 
tices [llMlq . These systems have been experiment ally 
realized by pumping bosons into high orbital bands |l7l - 
K^ . Bosons have been observed to develop phase coher- 
ence before they decay to the lowest band. In addition, in 
the artificial lattice systems of exciton-polariton in semi- 
conductor quantum wells, a d-wave condensation in the 
excited bands of bosons have also been observed (20| . 
These UBECs are beyond Feynman's "no-node" theorem 
0, im , which states that the ground-state wavefunctions 
of bosons are positive-definite under very general circum- 
stances. This theorem applies for the system of superfluid 
''He [2lj and many experiments of alkali bosons [l^]. It 
also implies that time-reversal (TR) symmetry cannot be 
spontaneously broken in usual BECs. However, UBECs 
escape from the "no-node" constraint. Their condensate 



wavefunctions are nodal, which are able to break TR 
symmetry spontaneously under certain conditions @. 

The recent UBECs realized in Hemmerich's group is 
an exciting progress [18| , where the time-of-flight (TOF) 
spectrum has revealed signatures of both the real and 
complex UBECs by tuning the anisotropy of the optical 
lattice. However, the TOF images can only provide the 
single-particle density distribution in momentum space, 
thus in the complex UBECs, the key information about 
the relative phase between the condensate components 
is lost during TOF. Without the phase information, the 
TOF images of the Px ± ipy BEC can be interpreted by 
other plausible scenarios such as the phase separation 
between two real condensates at different momenta or the 
incoherent mixing between them. It would be important 
to have the smoking gun evidence of the phase difference 
±-| between the two condensate components. 

In this paper, we propose a phase-sensitive detection 
scheme to identify the Px ± ipy symmetry of UBECs 
by measuring the relative phases of ib-^. This proposal 



is based on the scheme in Ref. [23|, which has been 
used to construct the off-diagonal correlation functions 
in momentum space for UBECs. By implementing a 
momentum-kick Raman pulse, we build up the connec- 
tion between bosons with different condensate momenta 
in the complex UBEC. As we will show below, the rela- 
tive phase information is uniquely tied to the off-diagonal 
correlations between the different condensate momenta, 
which can be measured in time-of-flight imaging through 
the impulsive Raman pulse. We note that a different 
scheme has been proposed recently by Kitagawa et al. 
for phase-sensitive detection of nontrivial pairing sym- 
metries in ultracold fermions based on the two-particle 
interferometry 24 1. 

Our scheme is connected with the experiment in Ref. 
[ist . where the bosons are pumped to the first excited 
band of a s-p hybridized system (hybridization between 
the s-orbital of the shallow sublattice A and the p-orbital 
of the deep sublattice B as illustrated in Fig. [Ija) . Two 
degenerate band minima (denoted as Ki,2 below) locate 
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FIG. 1: (a) The bipartite lattice structure (along the x- 
direction) in Hemmerich's experiment [l^, where bosons are 
loaded into the second band, (b) The first Brillouin zone (blue 
area) and the basis vectors of the reciprocal lattice 61,2- 



at the lattice momenta ^61^2 with bi^2 = (iTi": t^) the re- 
ciprocal lattice vectors (we set the lattice constant a = 1), 
as shown in FiglT](b). Bloch-wave states ^'1^2(1") at K12 
points are real- valued wavefunctions with nodal lines sim- 
ilar to standing waves, and thus are time-reversal invari- 
ant. The complex combination of = {"i^i ± i'^2)/V2 
only have nodal points from intersections of the nodal 
lines of ^1,2. With repulsive interaction, the complex 
condensates 'I'± with nodal points are favored since their 
spatial distributions of particle density are more uniform 
and extensive than other states, minimizing the interac- 
tion energy [lo| . 

Both of the wavefunctions ^'1^2(1") have odd parity with 
the p-wave symmetry. Rigorously speaking, the lattice 
configuration in Ref. [3] does not have 4-fold rotational 
symmetry, and thus ^'1^2 are not transformable to each 
other by the rotation of 90°. Nevertheless, for simplic- 
ity, we still denote these condensates with the Px± ipy 
symmetry. The TOF imaging in the experiment |18| has 
observed four peaks at {Ki^2, K[ 2 = -fi'1.2 — &1.2) with 
the same height, which implies that the condensate has 
equal weights of the components ^'i and ^£'2. However, 
the phase difference between the two components and 
^2, which is critical for verifying the novel ± ipy con- 
densation symmetry, is not clear. 




FIG. 2: Using two Raman pulses with different propagating 
directions to built up the momentum transfer between bosons 
at Ki and A'2. 



The spirit of our proposal for the phase-sensitive de- 
tection can be outlined of as follows: for the conden- 
sate \^r) = 



(C^, +e''Ciy^-\0), where 



the boson number in the condensate; C\^^ (C'lra) 
bosonic operator creating a boson at Ki (K2); is the 
relative phase between bosons in Ki and K2- If 6* 7^ 
and TT, the condensate exhibits a vortex-antivortex lat- 
tice structure. The staggered orbital angular momen- 
tum (0AM) density wave order parameter is defined 



L,(Q) = z(CT C 



. "^2 - CI^^Ck,) where Q ^ ~ K2 
Its magnitude reads as c\L z{Q)\^ c) = NcSinO, which 
is just the off-diagonal correlation. It reaches maximum 
for the Px ± ipy state (9 — i^^n). If we implement a 
Raman transition to transform the bosons in the original 
condensate into [2^] : 



K2 



V2 



{CK2-^CKJAl) 



we can measure the density difference of the new BECs 
Sn' = i^clC'lC'j,^ - C^^C^J^e) through the TOF 
imaging. Sn' exactly gives the desired off-diagonal corre- 
lation as 



Sn' ^t{^,\C],. C 



K2 



(2) 



Now we turn back to the Hemmerich's experiment and 
show how to implement the Raman transition. Similar 
to Ref. (2^, the Raman transition can be realized by 
two traveling-wave laser beams propagating along differ- 
ent directions (as plotted in Figl2]) with corresponding 
wavevector ki 2 and frequency LO12, which introduce an 
effective Raman Rabi frequency with a spatially varying 
phase n{v,t) = noe'^^^-''-^'^*+'t'\ where (5k = ki - ka, 
Slo=uji — UJ2, and </> is the relative phase between the two 
Raman beams, fio is expressed as flo = f^iilj/A, where 
A is the detuning, ^1(2) are the resonant Rabi frequencies 
for the individual transitions between the initial (finial) 
states and the intermediate state, and are proportional to 
the strength of the electric field of the corresponding Ra- 
man beams. This spatially dependent Raman transition 
builds up the connection between the condensation com- 
ponents in two degenerate points K12 = (±7r/2, 7r/2), 
which demands that Sk = K2 — Ki — (tt, 0) (Notice that 
^k and —Sk are connected by a reciprocal lattice vector 
and thus equivalent. Because of this feature, the Raman 



scheme here is simpler compared with the one in Ref. 23 1 
which needs to use transitions between two different hy- 
perfine levels) . The effective Hamiltonian for the Raman 
process is described by 



J drn{r, t)^"! {r)^{r) + h.c, 



(3) 



which, together with the original atomic Hamiltonian Hq 
in the optical lattice [18j, gives the full Hamiltonian of the 
system. 
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FIG. 3: The energy spectra along the lines between the high 
symmetry points with the experimental values of the param- 
eters given in Ref . . 

Different from Ref.^^, in tlie experiment 18] the op- 
tical lattice is too shallow to apply the tight binding ap- 
proximation. Instead, we expand the field operator in 
the Bloch representation as 

Vj/W-^CnkV'nkW, (4) 
nk 

where C„k is the bosonic operator annihilating a bo- 
son in the nth band with momentum k, tAnk(r) is the 
Bloch wavefunction, and the summation of k is over 
the first BZ. We choose the effective Rabi frequency f^o 
of the Raman pulse so that it is small compared with 
the band gap but large compared with the atomic hop- 
ping rate in the lattice {t ^ HHq ^ A 12). Under this 
condition, we can neglect the interband tunneling as 
well as the time-dependence of the wavepackets 5'i(r) 
and ^'2(r) during the Raman transition. For the typ- 
ical values of the experimental parameters, the energy 
band structure is shown in Fig. [31 The hopping rate 
is estimated by i « Q.ObEr ~ 27rfi, x 0.1 kHz and the 
smallest bandgap A12 w l.OSEr ~ 2nh x 2.2 kHz. If 
we choose hfta ~ 27r?i x 0.5 kHz, the corrections to 
the above approximation, estimated by t'^/ (?ii7o)^ and 
{hn.o)'^ / AI2, are pretty small. In our case, the Raman 
operation induces a transition between the complex and 
polar UBECs with the same kinetic energy but different 
interaction energy. hSuj should match the energy differ- 
ence between the initial (complex UBEC) and final (real 
UBEC) states of the Raman transition, which can be es- 
timated as 10~^Er ~ 2TTh x 2.1 Hz and much smaller 
than hflQ. Therefore, the phase accumulation induced 
by 6uj within the duration of Raman pulses 6t can be 
neglected {6u!6t ^ 7r/4), and the Raman Rabi frequency 
in Eq.([3l) can be considered as time- independent during 
the Raman transition. 

Under the above approximations, is simplified to 

i/H-e'*E^W^k+5k^^k + /i.c, (5) 

k 



FIG. 4: TOP imagings after the Raman transition with = 
for (a) the complex UBEC (^1 + ^^'2) and (b) incoherent 
mixing of the real UBECs (*i and ^2). 

which is constrained only to the relevant band. The sum- 
mation of k is over the first BZ, and the effective Raman- 
Rabi frequency fl(k) is k dependent and can be calcu- 
lated based on the eigenvectors obtained in the band- 
structure calculation. Before the Raman transition, the 
momentum distribution of the Px ± ipy condensate is 
sharply peaked at Ki and K2 with a small distribution 
width A, and within this small region J7(k) can be con- 
sidered as a constant, estimated by 17 ~ 0.98Jlo for the 
typical values of the experimental parameters. So dur- 
ing the Raman transition, the small wave-packets around 
Ki and K2 are transferred by the same formula without 
distortion. When we choose the duration of the Raman 
pulse as fl5t ~ 7r/4, the wave-packets around Ki and if 2 
with |(ik| < A are transferred by 

Cki+dk = (Cxi+dk - ie"^CK2+dk)/V2, 
Ck,+du - (C^K.+dk - *e-^^CK,+dk)/V2. (6) 

Right after this Raman transition, we turn off the trap- 
ping optical potential and perform the TOF imaging to 
measure the particle density distribution in momentum 
space around Ki(2). The density difference {Sn' (dk)) = 

('^K2-l-(ik~"A'i+dk) = (^ifa+dk^A'a+dk^C'ifi+dk'^A'i+dk) 

gives the off diagonal correlation for the original UBEC 

6n' (dk) = *(e^^C^^+,kCK,+,k - e-^^C%^^^CK,+dk) ■ 

(7) 

Notice that the experimental observable Sn' (dk =0), the 
height difference between the peaks in Ki and K2, is de- 
pendent on the phase difference of the two Raman pulses 
(j>. For the px ± ipy BEC, after the Raman pulse, we 
see that 5n' {dk —0) oc cos(f> from Eq.([7|). The oscilla- 
tion of Sn' with (p indicates coherence of the two Ra- 
man pulses, which is critical for phase-sensitive detec- 
tion. For = 0, Sn'{dk=0) = i{C]^^CK2 - C^kJ^k^) 
represents the order parameter of the orbital ordering of 
the original UBEC. With this phase-sensitive measure- 
ment, we can easily distinguish the complex condensate 
I*) « 2iVc/2vw (gAi+e''gL)'^'|Q) and other plausible 
scenarios, such as the phase separation or incoherent mix- 
ing between the polar UBECs l^-i) = -^(c]^J^<=|0) 

and 1 5-2) = 7^(CaJ^'=|0). In the conventional TOF 
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imagings, both of them exhibit four peaks with the same 
height thus can not be distinguished. Under this phase- 
sensitive TOF imaging, we find that for the ± ipy 
UBEC, (*kiK2+dkl*) = 0' which means that after Ra- 
man transition, the initial complex UBEC turns to the 
polar UBEC condensing only at Ki, and the predicted 
new TOF images are shown in Fig|4] (a) with only two 
peaks, in contrast with the four peaks that one expects 
to see for the incoherent mixing state between the con- 
densates 5'i and ^2 show in Fig. |4] (b). 
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FIG. 5: (a) The vortex-antivortex lattice pattern of the 
Px ± ipy BEG before the Raman transition; (b) the real-space 
current pattern during the Raman transition, solid (dashed) 
arrows denote original (reflected) currents. 

To get a better understanding of our results, we pro- 
vide a real space picture to illustrate the complex-real 
UBEC transition during the Raman process. Before 
the Raman transition, the Px ± ipy UBEC exhibits a 
vortex-antivortex lattice structure in sublattice B {p- 
orbital sites), as shown in Figl5l[a). The Raman beams 
introduce an extra potential with the form of Eq. ^ , as 
shown in FiglS] (b). Without loss of generality, we fo- 
cus on one site in sublattice B, initially the local wave- 
function within this site can be approximately considered 
as (/3(r) e'^[(^^(r) -|- i(py(r)], where (px(y)[Y) denote 
the p2,(j,)-orbital Wannier function. During the Raman 
transition, the initial current is reflected by the extra 
potential, and the local wavefunction turns to {p{t) = 
cos(fit)i^(r) -|-sin(rit)(p*(r), where </5*(r) is the TR coun- 
terpart of ip{v) carrying a current with an opposite direc- 
tion. Initially, </3(0) = ^piv) denotes the px + ipy state, at 
the momentum of t^Vl = tt/A, the reflected current hap- 
pens to cancel with the initial one, and it turns to a polar 
state with the real wavefunction </j(to) 'Pxi''^) — 'y'y(r), 
and the corresponding polar UBEC exhibits two peaks 
in the TOF spectrum, as shown in Fig|4] (a). 

In above analysis of the Raman transition, we have ne- 
glected the effect of interaction, which is responsible for 
the broadening of the TOF imaging peaks 2^, 2^ . Now 
we estimate the interaction effect by solving the time evo- 
lution from the Gross-Pitaevskii (GP) equations (the CP 
equation gives an adequate description of the interaction 
since the initial state of the system is a BEC). A simi- 
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FIG. 6: Time evolution of the condensation fraction at ^Pki 
with Raman frequency Q, = 2n x 0.5kHz, and interaction pa- 
rameter gop — (solid line) and O.GEr (dashed line), respec- 
tively 



lar problem has been addressed for the many-body Rabi 
oscillation in a two-component BEC As analyzed 

above, we neglect the deformation of ^'i(r) and ^'2(r) 
during the time evolution, and the dynamics of the sys- 
tem can be approximately considered as two-mode tran- 
sition, thus the wavefunction during the time evolution 
can be expressed as: *(t) = Ci(i)*i + C2(i)*2- The 
corresponding GP equation reads: 



= {Ho + n{r) + .gop|*(r0p}vl/(f). 



(8) 



dt 

where Hq = — -I- Vo(r) is the original optical lattice 
Hamiltonian, p is the average density, go is the s-wave 
scattering interaction parameter. In the experiment [ist . 
(7oP is estimated to be O.Gi^r- Notice that '!/'i(2) are eigen- 
functions of ifo- thus we get: 

^ ^7C2 + (25|Clp + 4g'|C2nCl + 25'C^C2^ 

.dC2{t) 



where g = g^p j (Pr\^i{r)\'^ = g^^p j cP'r\^2{r)\'^ , 
9' = 9oP J d'^r\^ii{r)\'^\'^2ir)\'^ .Using the initial condi- 
tion: Ci(0) = I/V2, 6*2(0) = i/V2 (^-Ki +i'^K2 state) 
and O = 27r X 0.5 kHz, we obtain the time evolution 
of condensation fraction of bosons at Ki wave-packet, 
and compare the result with the non-interacting case. 
As shown in Fig. [6l the interaction barely changes the 
Rabi oscillation within the duration of the Raman pulses 
in our case, which implies the single-particle Rabi os- 
cillation approximation we adopted above provides an 
accurate description for the many-body dynamics of the 
system during the Raman transition. 

To conclude, we propose a phase-sensitive detection 
scheme to identify the nontrivial symmetry of recently 
observed Px ± ipy orbital UBEC, where the connection 
between different condensate components is built up by 



5 



the momentum kick provided by impulsive Raman pulses. 
Our scheme can also be applied to the phase-sensitive de- 
tections for unconventional BECs with other symmetries, 
eg. the recently observed UBEC in /-orbital bands of the 
optical lattice jl9| . 
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